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1.  Introduction  Kettelle  (1962)  presents  an  algorithm  for  allocating 
redundancy  so  as  to  maximize  system  reliability  without  exceeding  a 
specified  cost  (or  equivalently,  achieve  a  specified  reliability  at 
minimum  cost) .  In  the  present  paper  we  develop  an  algorithm  to  solve 
the  more  general  problem  of  maximizing  system  reliability  without 
exceeding  any  of  several  linear  constraints. 

Specifically,  we  consider  a  system  consisting  of  k  "stages". 

The  system  functions  if  and  only  if  each  stage  functions.  Stage  i 

consists  of  n^  (to  be  determined)  units  of  type  i  in  parallel,  so 

that  stage  i  '  functions  if  and  only  if  at  least  one  of  the  n^  units 

of  type  i  function,  i  =  l,2,...,k.  Suppose  unit  i  has  a  "cost"  c^ 
th 

of  the  j  type,  i  =  l,...,k;j  =  l,...,r.  As  an  example,  the  first 
type  of  cost  might  be  weight,  the  second  volume,  and  the  third  money. 

A  linear  constraint  exists  on  each  cost  as  follows: 
k 

J^ifi  -  Cj,  j  =  l,2,...,r.  (1) 

Thus  in  the  example  the  total  weight  of  the  system  might  be  required  not 
to  exceed  a  specified  amount  c^,  the  total  volume  required  not  to  exceed 
Cg,  and  the  total  cost  in  dollars  required  not  to  exceed  c y 

A  unit  of  type  i  has  probability  p^^  of  functioning,  independently 
of  the  functioning  or  non-functioning  of  the  other  units  of  the  system. 
Thus  system  reliability  P(n) ,  where  n  =  (n^...^),  is  given  by 
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where  q^  =  1  -  p^.  Our  problem  is  to  choose  n  (a  vector  of  positive 
integers)  so  as  to  maximise  P(n)  in  (2),  subject  to  constraints  (1). 

2.  Domination  Let  c,(n)  =  Z  c.  .n.  represent  the  j  1  cost  of 

J  “*  4  — T  *  J  *** 

±~~±  i  2 

the  redundancy  allocation  n.  Then  n  is  said  to  dominate  n  if 

c^(n^)  <  Cj(n^),  3  =  while  P(n'*')  ^P(n^).  If,  in  addition, 

1  2 

at  least  one  inequality  is  strict,  then  n  is  said  to  dominate  n 
strictly.  A  sequence  S  of  redundancy  allocations  n*1,  h  =  1,2,..., 

satisfying  the  c ons  traints  (l)  , " Is . saTd’'  t6"be'X’aoMhatlhg“s¥qu®i5e"^ 

no  n1  is  strictly  dominated,  and  if  every  n  satisfying  the  constraints 
(l)  which  is  non  strictly  dominated  occurs  in  S. 

It  is  clear  that  to  solve  our  problem  we  need  only  consider  the 
members  of  the  dominating  sequence  S.  SpecificallyT  we  seek  that  alloca¬ 
tion  of  S  with  maximum  reliability  P(n)  among  the  members  of  S. 

Procedure  for  Two-Stage  System  First,  to  generate  the  dominating 

sequence  corresponding  to  a  system  consisting  only  of  stages  1  and  2,  we 

set  up  a  two-way  table  in  which  the  entry  in  row  n^,  column  n^  consists 

of  the  vector  (c^(n^,n2)  ,c2(n^,n2^ . . .  ,cr(n^,n2)  ,Q(n^,n2) ) ,  where 

^1  n2 

c^(n1,n2)  =  c1jn1  +  c^n^  j  =  l,...,r  and  QU^n^  =  1  -  (1-q.j.  )(l-q2  ))• 
This  is  the  vector  of  costs  and  unreliability  resulting  from  using  n^  \ 
units  of  type  1  and  n^  units  of  type  2.  Only  entries  satisfying  the 
constraints  (1)  are  included.  We  then  eliminate  from  the  table  any 
strictly  dominated  vector,  that  is  any  vector  all  of  whose  coordinates 
are  at  least  as  large  as  the  corresponding  coordinates  of  some  other 
vector  in  the  table,  with  strict  inequality  for  at  least  one  coordinate. 

The  remaining  undominated  allocations  constitute  a  dominating  sequence. 
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See  the  worked  example  of  Section  5  to  help  clarify  these  steps. 

Next  we  shall  show  that  we  can  construct  the  dominating  sequence 
corresponding  to  an  s  stage  system  from  the  dominating  sequence  corres¬ 
ponding  to  a  subsystem  of  s  -  1  stages.  Once  this  is  established  we 
will  then  be  able  to  construct  the  dominating  sequence  for  the  full  k 
stage  system  recursively;  i.e.,  first  for  a  subsystem  consisting  of  stages 
1  and  2,  next  combining  the  resulting  dominating  sequence  with  stage  3 
to  yield  the  dominating  sequence  for  stages  1,2,  ana  3  combined,  etc., 
until  the  dominating  sequence  for  the  full  k  stage  system  is  obtained. 

The  following  procedure  includes  the  Procedure  for  Two-Stage  System  as 
a  special  case. 

Procedure  for  s  Stage  System  (called  Procedure  for  short)  Set 

up  a  two-way  table  in  which  the  n  1  row  corresponds  to  n  units  of 

s  s 

type  s,  while  the  htu  column  corresponds  to  nh,  the  member  of 

the  dominating  sequence  for  the  first  s  -  1  stages.  The  entry  at  sow 
n  ,  column  h,  is  the  vector  (c  (nh  ,n  )  ,c_(nh  ,n  ),..., c(nn  ,n  )  ,Q(nl4,n  ) ) , 

S  i.*-  <srs  —  s 

the  vector  of  costs  and  unreliability  resulting  from  using  the  vector 

(n^,n  ).  Note  that  c.(nk,n  )  =  c.(n^)  +  c  ,n  ,  j  =  l,...,r,  while 
-  s  J  ~  3  n  J  “  sj  3 

Q(n  ,n  )  =  1  -  (1  -  Q(n  ))(1  -  q  3) .  Only  entries  satisfying  the  con- 
straints  are  included.  We  eliminate  from  the  table  any  strictly  dominated 
vector  (strictly  dominated  by  some  other  vector  in  the  table) .  The 
remaining  entries  constitute  the  dominating  sequence  for  stages  l,2,...,s, 
as  we  prove  in 

Theorem  1  The  vectors  that  remain  strictly  undominated  in  the 
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two-way  table  generated  in  the  Procedure  constitute  a  dominating 
sequence  foi  the  s -stage  system. 

Proof;  We  need  to  prove  (a)  the  allocations  obtained  following 
the  Procedure  include  all  strictly  undominated  allocations,  (b)  every 
allocation  obtained  using  the  Procedure  is  strictly  undominated. 

We  prove  (a)  inductively.  First  note  that  for  a  single  stage 

sys  tern all  “ aTToc a’€ibns^re"'sT^Tc^Iy^n3MInatg3d;  - 

the  allocations  obtained  by  the  Procedure  for  a  j  stage  system,  where 
j  =  l,2,...,s  -  1,  include  all  strictly  undominated  allocations  satis¬ 
fying  (1).  Consider  any  allocation  n  =  (n. , . . .  ,n  )  satisfying  (l). 

“X  s 

Then  by  inductive  hypothesis  (nn,...,n  ,)  is  dominated  by  some 

x  s  —x 

strictly  undominated  allocation  (n#,..,,n*  ,)  obtained  by  the  Procedure. 

Thus  by  definition  Q(n1, . . .  ,ns_1)  £  Q(n£, . . .  ,n*_1)  ,c  j(n1» . . .  jn^) 

>  c.(n#,...,n#  . ) ,  j  =  1, . . . ,r.  It  follows  that  Q(n)  =  1  -  P(n, ,...,n  ,)P(n 

J  X  S  —X  X  S  —X 

>  1  -  P(n*,...,n#  , )P(n*)  =  Q(n#)  ,  where  n#  =  n  ,  and  that  c.(n) 

““  XS-XS—  S3  J  ”* 

=  c<.(n1,...,ns_1)  +  Cj(ng)  >  c^(n* , . . .  .n*^)  +  c..(n#)  =Cj(n*),  j  =  l,...,r, 
so  that  n  is  dominated  by  n*.  On  the  other  hand,  n*  being  an  entry 
in  the  two-way  table  generated  by  the  Procedure  is  itself  dominated  by 
an  allocation  obtained  following  the  Procedure.  Thus  we  have  proved 
that  every  allocation  satisfying  (l)  is  dominated  by  some  allocation 
generated  following  the  Procedure.  Hence  the  inductive  proof  of  (a)  is 
completed. 

To  prove  (b) ,  suppose  n°  is  an  allocation  obtained  using  the 
Procedure.  If  n°  is  strictly  dominated  by  any  allocation  satisfying 
(l)  it  must  also  be  strictly  dominated  by  some  undominated  allocation 
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satisfying  (l).  But  we  have  just  proved  that  all  undominated  allocations 
satisfying  (l)  are  obtained  by  the  Procedure.  Thus  n°  is  strictly 
dominated  by  say  n^  also  obtained  by  the  Procedure.  This  is  a  contra¬ 
diction  since  no  allocation  obtained  under  the  Procedure  can  dominate 
any  other  allocation  obtained  under  the  procedure. 

3*  Approximations  In  applications  of  the  Procedure  we  may  generally 

-apply- the.  following  . approximation. _ Instead  of  using _ 

nn  n.  n..  n_  n..  n_ 

Q(n1,n2)  =  1  -  (1  -  q1~)(l  -  qg  )  =  q^1  +  q,  -  qx  q2  , 

we  disregard  the  product  term  and  use 

nl  n2 
Q(n1}n2)  -  q^  +  q2  . 

Similarly,  for  an  s  stage  system,  we  approximate  (where  n  =  (n.,,...,n  ..) 

a-  S  —J. 

Q(n,ns)  =  Q(n)  +  qg .  (3) 

Kettelle  (1962)  shows  that  the  use  cf  this  approximation  for  the  case 
r  =  1  results  in  an  error  in  P,  the  system  reliability  achieved,  of 
at  most  Q  (where  P  +  Q  =  l) .  For  the  present  case  of  r  >  1,  the 
proof  goes  through  just  as  in  Kettelle  (1962).  We  do  not  repeat  the 
details . 

In  any  applications  of  the  Procedure  we  will  use  approximation 
( 3)  throughout . 

Another  approximation  that  may  reduce  the  length  of  dominating 
sequences  is  the  following.  In  comparing  a  pair  of  entries  in  the  table 
developed  following  the  Procedure  we  may  introduce  a  tolerance  factor 
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e  for  the  j  n  type  of  cost  and/or  a  tolerance  factor  e  for  unrelia- 

XU 

bility.  If  two  entries  in  the  table  differ  by  or  less  in  the  j 

th 

type  of  cost,  they  are  considered  alike  as  far  as  the  j  type  of  cost 
is  concerned;  similarly  if  they  differ  by  or  less  in  unreliability. 
The  result  is  that  domination  becomes  more  likely  so  that  the  lengths 
of  dominating  sequences  are  reduced.  Problems  that  might  be  otherwise 
unsolvable  because  of  excessively  long  dominating  sequences  can  some¬ 
times  be  solved  (approximately)  by  introducing  one  or  more~toTerince 
factors.  The  most  conservative  procedure  is  first  to  try  to  solve  the 
original  problem  without  tolerance  factors.  Then  if  the  dominating 
sequences  are  too  long  to  permit  a  solution,  introduce  a  small  tolerance 


factor  on  the  unreliability.  If  the  dominating  sequences  are  still 

too  long  either  increase  e  or  introduce  additional  tolerance  factors 

q 

Continue  increasing  or  adding  tolerance  factors  until  a  solution 
is  attained. 


4.  Starting  Values  for  the  n^,.  As  we  shall  see  later  the  lengths 
of  the  dominating  sequences  determine  the  limiting  size  of  the  problems 
that  can  be  handled  on  a  computer  and  the  time  required  to  compute  solu¬ 
tions.  It  is  therefore  of  utmost  importance  that  the  lengths  of  the 
dominating  sequences  be  kept  as  small  as  possible.  One  way  to  reduce 
their  size  is  to  use  as  large  a  starting  value  for  each  n^  as  is 
conveniently  possible. 

We  now  describe  a  method  for  finding  such  large  starting  values. 

(1)  Add  one  unit  of  each  component  type  in  succession  until  finally  a 
constraint  will  be  violated  upon  the  next  addition.  (2)  Compute  the 
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reliability  P  of  the  resulting  system.  (3)  Determine  n^,  the  minimum 
number  n  of  unite  of  type  i  required  to  achieve  a  reliability  of  P 
or  greater,  from 

P  <  1  -  qj*  (4) 

Then  it  is  clear  that  the  solution  to  the  allocation  problem  will  require 

a  value  of  n^  at  least  as  large  as  n°.  (4)  Thus  the  starting  value 

. -  - -  - .  . . .  & 

of  n^  may  be  taken  as  .  . .  ~ 

To  illustrate  the  value  of  this  starting  procedure,  note  that  in 
on?  problem  involving  10  stages  and  3  constraints  using  this  procedure 
reduced  the  length  of  the  dominating  sequence  from  the  starting  point 
until  a  constraint  was  violated  from  334  to  62.  In  a  second  problem 
involving  20  stages  and  3  constraints,  without  this  procedure  the  computa¬ 
tion  had  to  be  halted  at  the  10^  stage  with  the  length  of  the  dominating 
sequence  559,  while  the  use  of  the  procedure  led  to  a  solution  with  the 
final  dominating  sequence  of  length  69. 

An  alternate  method  for  generating  starting  values  for  the  n^  is 
to  use  tolerance  factors  as  described  in  Section  3  to  obtain  an  approxi- 
matr  solution.  After  the  approximate  solution  is  obtained,  use  steps 
(2), (3),  and  (4)  above. 

5.  Example  We  illustrate  the  application  of  the  algorithm  with  the 
following  worked  example,  based  on  one  that  appeared  in  Kettelle,  1962. 
Consider  a  four-stage  system  with  tinlt  costs  and  unreliabilities  as 
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follows i 


Stage,  i 

1 

2 

3 

4 

Constraint 

Cost,  c^ 

1.2 

2.3 

3.4 

4.5 

47.0 

Weight,  ci2 

1 

1 

1 

1 

20 

Unreliability, 

«i 

.2 

.3 

.25 

.15 

wish  to  choose 

(n^n^n^n^) 

so  as 

to  maximize  system  reliability 

(2)  explicitly  given  by 

P(n;L,n2,n3,n4)  =  (1  -  ^(l  -  ^(l  -  ^^(l  -  .15%,  (5) 

subject  to  constraints  (l)  explicitly  given  by 

1.2nx  +  2.3n2  +  3.4n.j  +  4.5n4  <  47.0, 

(6) 

nl  +  n2  +  n3  +  n4  £  ^0. 

First  we  shall  obtain  starting  values  for  the  n^  following  the 
procedure  of  Section  4*  Starting  with  n^  =  l,n2  =  l,n3  =  l,n4  =  1,  we 
add  one  component  at  a  time  until  adding  an  additional  component  would 
violate  a  constraint.  As  shown  in  Table  1  we  arrive  at  a  system  composi¬ 
tion  (5, 4, 4, 4)  with  system  reliability  .9872.  Solving  for  n°,  the 
minimum  n  satisfying 

1  -  .2n  £  .9872, 

we  obtain  n^  =  3«  Similarly  n2  =  4,^  =  4,n°  =3,  as  shown  in  Table  2. 
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Table  1  -  Computation  to  Find  an  Attainable  Reliability 
before  Violating  a  Constraint 


Stage 

2 

Stage 

? 

Stage 

4 

Cost, 

J=1 

Weight 

-1=2 

1 

1 

l 

11.4 

4 

1 

1 

l 

12.6 

5 

2 

1 

l 

14.9 

6 

Constraint  1  =  47.0 

2 

2 

l 

18.3 

2 

2 

2 

22.8 

Constraint  2  =  20 

2 

2 

2 

24.0 

3 

2 

2 

26.3 

10 

3 

3 

2 

29.7 

11 

3 . 

...JL  _ _ 

. 3 _ 

JA*2 . 

12 . 

3  1 

3 

3 

35.4 

13 

4 

3 

3 

37.7 

14 

4 

3 

41.1 

15 

4 

4 

45.6 

16 

4 

4 

46.8 

17 

<-  Attainable  Reliability  =.9872 

5 

4 

49.1 

18 

■ 

Table  2  -  Computation  of  Minimum  No.  of  Each  Stage  Required 
to  Achieve  Attainable  Reliability  of  .9872  for 
That  Stage  Alone 
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Next,  following  the  Procedure  of  Section  2,  we  obtain  the  dominating 
sequence  for  Stages  1  and  2  combined,  as  shown  in  detail  in  Table  3. 


Table  3  -  Dominating  Sequence  for  Stages  1  and  2 


Note  that  we  have  listed  the  costs  and  unreliability  starting  with 
n°  =  3  and  adding  one  unit  at  a  time  for  Stage  1  across  the  top. 
Starting  with  n°  =  4  and  adding  one  unit  at  a  time  we  list  the  costs 
and  unreliability  for  Stage  2  down  the  side.  We  obtain  entries  in  the 
body  of  the  table  by  adding  the  respective  costs  and  unreliabilities} 
only  entries  satisfying  the  constraints  are  retained.  Proceeding 
systematically,  comparing  pairs  of  entries,  we  eliminate  all  strictly 


| 


11 


dominated  entries.  Thus  the  entry  in  row  1,  column  4  is  eliminated 
since  it  is  daminatsd  by  the  entry  in  row  2,  column  2.  (Note  that 
each  figure  in  the  latter  position  is  less  thsn  the  corresponding 
figure  in  the  former  position.)  Similarly  the  entries  shown  with  a 
line  through  them  are  strictly  dominated.  The  remaining  entries  are 
not  strictly  dominated.  Note  that  only  a  portion  of  the  complete  table 
is  shown;  actually  entries  continue  to  be  made  until  a  constraint  is 
violated . 

In  Table  4  the  dominating  sequence  is  obtained  for  Stages  1,  2, 
and  3  combined.  Across  the  top  of  Table  4  are  listed  the  entiles  of 
the  dominating  sequence  of  Stages  1  and  2  obtained  in  Table  3>  Down 
the  side  are  listed  entries  corresponding  to  adding  one  unit  of  type  3 
at  a  time  starting  with  n°  =  4.  Chly  entries  satisfying  the  constraints 
are  listed  in  the  body  of  the  table.  Again  proceeding  systematically 
comparing  pairs  of  entries,  we  strike  out  all  dominated  entries.  Thus 
the  entry  in  row  2,  column  2  is  eliminated  since  it  is  dominated  by 
the  entry  in  row  1,  column  4>  The  remaining  entries  constitute  an 
undominated  sequence. 

Finally  we  form  Table  5  to  yield  the  dominating  sequence  for 
Stages  1,  2,  3,  and  4  combined.  Proceeding  as  before  we  obtain  the 
dominating  sequence  for  the  full  system.  The  solution  to  our  problem 
is  the  entry  in  the  table  with  lowest  unreliability,  namely  the  entry 
with  costs  c1  =  46.8,  c2  =  17,  and  unreliability  .00840.  To  obtain 
the  corresponding  system  composition  we  must  trace  back  through  Tables  3, 


.  asa  aaa 
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Table  5  -  Dominating  Sequence  for  Stages  1,  2,  3  and  4 


4,  and  3.  Mote  that  from  Table  5  the  optimal  n^  =  3,  while  from 
Table  4,  the  optimal  =  5.  From  Table  3,  the  optimal  n^  =  4, 

=  5.  Actually  in  the  machine  program  (see  Section  6)  the  composition 
of  the  system  is  retained  at  each  stage  so  that  no  retracing  is  necessary. 
Thus  our  solution  consists  of  a  composition  of  n^  =  4>  n^  =  5,  n^  =  5, 

n^=  3,  with  associated  reliability  (exact  value)  from  (5)  of 

13(4,5,5,3)  =  (1  -  .24)(1  -  .35)(1  -  .255)(1  -  .153)  =  .9916. 

2 

Note  that  the  error  using  approximation  (3)  throughout  is  <  .0064  -.000071 

6.  Computer  Program  The  procedure  described  in  Section  2  above  was 

programmed  for  the  I3M  7090  Data  Processing  System.  The  program  was 
written  with  the  capacity  to  handle  one,  two,  or  three  cost  constraints, 
a  maximum  of  sixty-four  stages,  a  maximum  of  ten  units  for  each  stage, 
and  a  maximum  of  1024  entries  in  the  dominating  sequence  at  any  combining 
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step.  The  essential  features  of  the  program  are  presented  in  the  flow 
chart  of  Table  6.  The  following  notation  is  used  in  the  flow  chart* 


Input  Quantities:  k  -  no.  of  stages 


r  -  no.  of  cost  constraints 


°U 

pi 


-  value  of  the  j  “  cost  constraint 

til  til 

-  j  1  cost  of  one  unit  of  the  i  stage 

til 

-  reliability  of  one  unit  of  the  l  stage 

-  tolerance  for- unreliability  ~ 

-  tolerance  for  cost 


xu 

Computed  Quantities:  -  unreliability  of  one  unit  of  the  il  stage 

til 

-  no.  of  units  of  l  1  stage  to  start  algorithm 
PQ  -  trial  system  reliability 

til 

-  no.  of  row  headings  for  the  i  stage 

NCih  -  no.  of  units  of  the  ith  stage  at  the  hth 
column  heading 

M  -  no.  of  column  headings 


I  -  index  for  the  stage  currently  being  combined 
into  the  system 

QCH,  -  unreliability  of  the  system  represented  at 
the  hth  column  heading 

XU 

QRH^  -  unreliability  of  the  no.  of  units  of  the  I 
stage  represented  at  the  h^*1  row  heading 

tU  AL 

CCH.,  -  j  cost  of  the  system  represented  at  the  h 
column  heading 

XU  XU 

CRH..  -  j  n  cost  of  the  no.  of  units  of  the  I  stage 
J  represented  at  the  hth  row  heading 

XU 

TE. .  -  indicator  for  the  table  entry  at  the  i hn  row 
J  and  column  (0  indicates  entry  is  in  the 
dominating  sequence j  1  indicates  entry  is  not 
in  the  dominating  sequence) 


TNClh  -  temporary  storage  for  the  no.  of  units  of  the 
ith  stage  at  the  h^  column  heading 


The  circled  numbers  appearing  by  the  boxes  of  the  flow  chart 
refer  to  the  following  explanatory  notes* 

1.  The  input  quantities  required  are  entered  into  the  computer 
through  the  appropriate  input  device. 

2.  The  unreliability  of  one  unit  of  each  stage  is  computed.  The 
si  are  initialized  to  one. 

3.  This  loop  adds  one  unit  of  each  stage  successively  until  one 

of  the  cost  constraints  is  exceedecL  ~~  “ 

4.  The  last  s^  to  which  a  unit  was  added  is  reduced  by  one  so 

til 

that  the  system  represented  by  units  of  the  i  stage, 
i  =  l,...,k,  is  one  which  violates  no  constraints. 

5.  The  trial  reliability,  Pq,  of  this  system  is  computed. 

6.  The  s^  are  initialized  to  zero. 

7.  This  loop  computes  the  number  of  units,  s^,  of  each  stage 
required  so  as  to  make  the  reliability  of  that  stage  at  least 
equal  to  the  trial  reliability,  P^.  The  resulting  s^  are 
the  values  used  to  start  the  algorithm. 

8.  The  number  of  row  headings  for  the  i^  stage  is  computed  to 
be  the  largest  integer,  H^,  such  that  £  the  maximum 
number  of  units  permitted  for  each  stage,  and  n  •  c^.  £  c. 
for  j  =  1,. .  ,r. 

9.  The  units  count  for  the  first  stage,  NC^>  is  set  at  each 
column  heading,  h. 

10.  The  index,  I,  of  the  stage  to  be  added  to  the  system  is  set 
to  two.  The  number  of  column  headings,  M,  is  set. 

11.  The  unreliability  of  the  system  represented  at  each  column 
heading  is  computed. 


12.  The  costs  of  the  system  represented  at  each  column  heading 
are  computed. 

13.  The  unreliability  of  the  units  represented  at  each  row  heading 
is  computed. 

14.  The  costs  of  the  units  represented  at  each  row  heading  are 
computed. 

15.  This  loop  computes  the  costs  that  occur  at  each  entry  in  the 

. . tahla-iy-.addinff-.tha  nn sta  at.  the  respective  rnu  and  column _ 

headings.  If  any  of  the  costs  exceed  a  cost  constraint  the 
entry  is  marked  as  not  being  in  the  dominating  sequence, 
otherwise  the  entry  is  marked  as  being  in  the  dominating 
sequence. 

16.  This  loop  steps  the  indices  from  table  entry  to  table  entry 
and  checks  the  table  entry  indicator  to  see  if  the  entry  is 
in  the  dominating  sequence.  If  it  is  not, then  the  indices 
are  stepped  to  the  next  table  entry.  If  it  is,  then  the  next 
loop  is  executed. 

17.  This  loop  compares  the  table  entry  just  chosen  with  every 
other  entry  in  the  dominating  sequences  and  marks  every  entry 
which  it  dominates  appropriately.  It  then  goes  back  to  con¬ 
tinue  the  loop  described  in  16  above. 

18.  When  the  loop  described  in  16  above  has  been  exhausted,  the 
table  entries  marked  with  a  zero  are  the  true  dominating 
sequence.  The  units  counts  for  the  column  headings  are  now 
moved  to  temporary  storage  locations  so  that  the  column  headings 
may  be  updated.  . 


19 .  This  loop  update*  the  unite  counts  for  the  column  headings 

20.  If  all  stages  have  not  been  combined  into  the  system,  then 
the  index  I  is  stepped  for  the  next  stage,  the  number  of 


column  headings  is  set,  and  the  program  is  repeated  from 
point  A  . 

21.  When  all  stagea  have  been  combined  into  the  system,  this 
loop  selects  the  set  of  units  counts  having  the  maximum 

reUaBilitjrT" Tfregg~"unl'ta  counts  are  the  elements  of  the - 

vector  n  and  represent  the  optimum  redundant  system. 

22.  The  vector  n  is  printed  through  the  appropriate  output 
device . 

7.  Problem  Solving  Experience 

The  purpose  of  this  section  is  to  present  some  of  the  experiences 
which  have  resulted  from  attempts  to  solve  problems  of  varying  size 
and  complexity  using  the  computer  program  described  in  the  previous 
section.  This  discussion  is  included  because  it  is  often  difficult  to 
predict  whether  the  solution  to  a  given  problem  is  practical  without 
prior  experience  on  a  similar  problem.  Even  though  the  solution  to  a 
problem  may  appear  to  be  practiced  in  the  sense  that  the  basic  quantities 
do  not  exceed  the  program  limitations,  it  may  turn  out  to  be  impractical 
because  the  dominating  sequences  developed  exceed  computer  capacity.  It 
is  also  possible  that  the  computer  time  required  may  be  excessive.  The 
experience  recorded  here  may  provide  some  guideposts  in  estimating  the 
practicality  of  attempting  to  solve  future  problems. 


for  h=l,...*Mjj=l, 


TKC 


Table  6  -  Flow  Chart  of  Computer  Program  (con 
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A  series  of  problems  vss  devised  to  test  the  program  systematically. 
Table  7  presents  the  basic  data  for  the  problems  and  indicates  the  nature 
of  the  results.  One  additional  problem  involving  twenty-five  stages  was 
attempted  but  was  not  successfully  completed  because  the  dominating 
sequences  became  so  large  that  the  computer  time  required  vma  excessive. 
This  does  not  necessarily  mean  that  all  problems  involving  more  than 
twenty  stages  are  impractical.  The  only  positive  method  for  determining 
whether  it  is  practical  to  solve  a  given  problem  is  to  attempt  to  find 
the  solution. 


Table  7  -  Problems  and  Results 


No.  of 
Stages 

No.  of 
Constraints 

Upper  Limit 
Set  on  No. 
of  Units  of 
Each  Stage 

Tolerance 
Factor,  e 

Average 
Length  of 
Dorn.  Seq. 

Maximum 
Length  of 
Oom.  Seq. 

Error  in  P 

10 

3 

6 

0 

62 

113 

0 

10 

3 

6 

10“V 

35 

62 

0 

10 

3 

6 

10”6 

21 

47 

0 

10 

3 

6 

10"5 

9 

.  19 

0 

10 

3 

6 

O 
•  H 

2 

3 

17  x  10“5 

20 

2 

5 

0 

196 

341 

0 

20 

2 

5 

10”  7 

174 

155 

0 

20 

2 

5 

10”6 

27 

48 

63  x  10"8 

20 

2 

5 

10“5 

15 

21 

10"5 

20 

3 

5 

0 

124 

214 

0 

20 

3 

5 

10~v 

98 

150 

0 

20 

3 

5 

10"6 

64 

112 

0 

20 

3 

5 

10"5 

31 

58 

78  x  10”7 
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8.  Otjher  Problem  Solvable  by  the  Same  Method  The  method  developed 
above  applies  to  a  number  of  problems  other  than  the  redundancy  alloca¬ 
tion  problem. 

(l)  Spare  Parts  Kit  (See  Proschan,  I960.)  A  system  is  required 
to  operate  during  [0,t].  When  a  component  fails,  it  is  instantly 
replaced  by  a  spare,  if  one  is  available.  The  components  considered 
operate  independently  and  are  essential  to  continued  system  operation, 
so  that  a  shortage  of  any  component  results  in  system  shutdown.  Only 
the  spares  originally  provided  may  be  used  for  replacement. 

Let  P^(n)  be  the  probability  that  n  or  fewer  spares  of  type  i 
are  required  (i.e.,  n  or  fewer  failures  of  type  i  occur  during  [0,t]), 
i  =  l,...,k.  Then  the  probability  P(n)  of  system  survival  during  [0,t] 
if  a  spares  kit  of  composition  n  =  (n^,... 
k 

p(s)  =  IT  ^(n.) . 

i=i  1  1 

The  problem  is  to  choose  n,  a  vector  of  positive  integers,  so  as  to 
maximize  P(n)  subject  to  linear  constraints  (l). 


,11^)  is  provided  is  given  by 


Note  that  in  the  typical  application  the  failure  distribution  for 


component  type  i  is  often  taken  to  be  exponential,  1 

such  cases  P^(n)  is  the  Poisson  distribution 

A.t  (\.t)n 

P4(n)  =  e  (l+\it+***+ — — ■). 


(2)  Optimum  Spares  Kit  when  Repair  Is  Allowed  As  above,  the 
system  is  required  to  operate  during  [0,t].  When  a  component  fails 
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it  is  replaced  by  a  spare,  if  available.  Repair  is  begun  immediately 
on  the  failed  component.  There  are  units  of  the  i  type  simul¬ 
taneously  in  operation  in  the  system  with  n^  spares  available  for 

replacement,  i  =  l,2,...,k.  The  failure  distribution  of  a  component 

At 

of  type  i  is  exponential,  1  -  e  ,  while  the  repair  distribution 
is  arbitrary  with  mean  with  w^X^pi^  <  1.  All  failure  and  repair 
times  are  independently  distributed. 


The  system  is  considered  to  have  failed  if  for  any  component  type 
no  spare  is  available  to  replace  a  failure;  i.e.,  if  say  for  type  i, 
a  failure  of  one  of  the  w^  operating  components  occurs  while  all  n^ 
spares  are  simultaneously  under  repair.  It  can  be  shown  (Karush,  1957) 
that  under  these  assumptions,  the  steady  state  probability  P(n)  that 
the  system  will  be  "available"  (i.e.,  not  shut  down  due  to  shortage) 
is  given  by 

k 

P(n)  =  TT  (nj) , 

i=l  1  1 


where 


Pi(ni)  = 


ni 

2 

h=0 


(Wi^l^i) 


hi 


h  ,  n^+1 

~/i  o 


(wixi^i) 

hi 


(8) 


As  before,  the  problem  is  to  choose  n  a  vector  of  positive 
integers  so  as  to  maximize  P(n)  subject  to  linear  constraints  (1) . 
Karush,  1957,  shows  how  to  solve  (approximately)  the  problem  when  a 
single  constraint  is  present  (r  =  l). 


The  algorithm  presented  above  applies  to  the  solution  of  both 
problems  (l)  and  (2).  The  only  change  is  to  use  the  P^(n^)  appropriate 
for  the  particular  problem,  in  carrying  out  the  Procedure  of  Section  2. 
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